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WEIGHTED TV MINIMIZATION AND APPLICATIONS TO VORTEX 

DENSITY MODELS 


P. ATHAVALE, R. L. JERRARD, M. NOVAGA, AND G. ORLANDI 

Abstract. Motivated in part by models arising from mathematical descriptions of 
Bose-Einstein condensation, we consider total variation minimization problems in which 
the total variation is weighted by a function that may degenerate near the domain 
boundary, and the fidelity term contains a weight that may be both degenerate and 
singular. We develop a general theory for a class of such problems, with special attention 
to the examples arising from physical models. 


1. Introduction 

Consider the weighted Total Variation functional 

Jiu)= [ f{x)\Du\:= sup I / u-div(/(x)C)dx : C e 1^1 < > (1-1) 

Jn IJn J 

defined for maps u G C M” open bounded. We always assume that the weight 

/ is positive and continuous in Q and possibly vanishing on some portion of dfl. 

We are interested in the analysis of functionals of the type 

h(u) = J(u) + ^ [ \u-uofg{x)dx (1.2) 

2 Jn 

with A > 0, g{x) > 0 a measurable weight function, and uq G := g{x)dx). 

Notice that, when f{x) = g{x) = 1, I\ corresponds to the well-known Rudin-Osher- 
Fatemi Total Variation based denoising model [25]. Anisotropic versions of TV functionals 
and applications to active contour and edge detection have been studied in [15]. Related 
homogeneous functionals for the description of landsliding have been also studied [20, 
13]. This model can be interpreted by convex duality as a kind of non-local vectorial 
generalization of the classical obstacle problem. Other vector-valued generalizations can 
be considered, such as J{u) = ||T • Du||, T a tensor (see for example [12] for the case 
T ■ P = trace (P), related to a Hele-Shaw model, corresponding to J(n) = ||divu||), and 
connections with game theory (tug-of-war type games) of the dual variational problem 
are possible. 

Due to the variety of applications of this model we are interested in a systematic study 
of properties of minimizers, in particular criteria for the minimizer being nontrivial, the 
characterization of the coincidence set, and conditions yielding the presence of flat zones 
and/or free boundaries. 

We are particularly interested in functionals of the type of /y that arise as reduced 
models in superconductivity and superfluidity (namely Bose-Einstein condensation) in 
axisymmetric 3-D domains (see e.g. [9, 10]), suited to describe the distribution of vorticity 
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that takes place in certain asymptotic regimes. Indeed, it is proved in [9, 10] that in certain 
scaling limits, a Bose-Einsten condensate may be described by a velocity field p : 
where H is a bounded, open subset of M^, that minimizes a functional of the form 

£{y) . / Pj,p 

Jv 

Here 

• Ppp is a nonnegative integrable function known as the “Thomas-Fermi density”, 
typically smooth in the interior of its support, corresponding to the limiting con¬ 
densate density in the asymptotic regime under consideration. 

• P := {x G : Pj,p{x) > 0 } is the limiting region occupied by the condensate 

• <I> is a given vector field that models forcing terms. 

• V X u is interpreted as the vorticity. In general it is only a measure, and | V x u| 
must be understood as a total variation. 

The model case occurs when Ppp{x) = (constant—quadratic)^ and <h(x) = a(—X 2 , xi, 0) 
for some a G M, corresponding to a quadratic trapping potential and forcing via rota¬ 
tion around the X 3 axis. In [9] it is noted that for this <h, if p is rotationally symmetric 
about the xs-axis (that is, p^p has the form Ppp{x) = Ppp{r,x^) for r = \/x\ + x|), then 
the minimizer has the form u* = rc*(r, X 3 )(^§^, ^, 0 ), for Wmin minimizing the reduced 
energy 

j Ppp(^\Dw\ + — — ^ ^ dr dz, (1.3) 

where 

n := {{r,z) G (0,oo) x M : Ppp{r, z) > 0}. 

Here, vortex lines are level sets of tc*, and “flat spots” - sets of positive measure on which 
is constant — represent regions in which the flow is irrotational. If H is bounded, then 
this can be written as a functional of the form I\, with coefficient functions / = Pj,p,g = 
“^Ptf I ^' 

All our results are formulated with this application in mind. In particular, we will 
always assume that 

/ G C(fl) nL“(fI), />0inH, (1.4) 

geLZm, g>0mn. (1.5) 

Note that we do not assume that /, g are bounded away from zero, or that g is integrable. 
(We require g G L^(H) only once, in Lemma 4.2, and it will turn out that other mitigating 
factors make the results of this lemma available for functionals of the form (1.3), even 
when g = Pj.p/r is not integrable.) The second half of this paper focuses on the specific 
example (1.3), for which the corresponding f,g described above satisfy (1.4), (1.5). 

The book [1] is a good general reference for mathematical issues relating to vortices 
in Bose-Einstein condensates. Most prior work on 3d condensates has dealt with an 
asymptotic regime in which the number of vortex filaments is bounded. As proposed in 
[3], then justified first in [19] and in greater generality in [22], in this limit the geometry 
of vortex filaments is governed by a limiting functional that is closely related to what one 


|V X u| -|- -|x — 
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obtains by considering in E above only the terms that are homogeneous of degree 1 with 
respect to v. Studies of this limiting functional are carried out in [2] and [23]. 

Plan of the paper: in Section 2 we revisit the theory of weighted BV functions and 
weighted TV functionals so as to include possibly degenerate weights, and we study 
general properties of minimizers such as existence, uniqueness and the Euler-Lagrange 
equation. 

In Section 3 we investigate properties of level sets of minimizers. Our main result in 
this section is Proposition 3.1, where we show a uniform lower bound on the measure 
of superlevel sets. As a consequence, for bounded minimizers, the maximal level set is 
flat and has positive measure. The proof is based on a suitable weighted isoperimetric 
inequality which may be of independent interest. 

In Section 4, under an additional regularity assumption on the weights, and using some 
ideas from [23], it is shown that minimizers are in fact of class BV (Lemma 4.1), and in 
dimension n = 2 level sets meet orthogonally the boundary (Lemma 4.2). 

In Section 5 we apply the previous theory to Bose-Einstein condensation models (for 
which g ^ Our main result in this section is Theorem 5.1, where we analytically 

derive properties of minimizers and their level sets, which represent vortex lines in the 
considered models. In particular, our results imply that in the limit in which condensates 
are described by the functional from (1.3), a condensate in the ground state always 

possesses a zone that vortx filaments do not penetrate, and in which the corresponding 
flow is thus irrotational. 

In Section 6 we show some numerical simulations, and eventually in the appendix we 
prove the weighted isoperimetric inequality used in Proposition 3.1. 

2. Weighted TV regularization - general results for rough weights 

2.1. Preliminaries on BV functions. The natural space for minimization of Ix is the 
weighted space 

H := L?g{EL) := {r : n —)■ M measurable : ||m|| < oo} 
with the norm and inner product 

guv dx. (2-1) 

Note that the weight g is not explicitly indicated in the notation for || • || or (•, •). With 
minimization problems in mind, we always work in the space H, even when discussing 
the pure weighted total variation J(-), without the quadratic term. Thus, all definitions 
that follow (subgradient, Legendre-Fenchel transform, etc.) are understood with respect 
to the H = Lg inner product. 

We also explicitly point out that, in this spirit, the domain of J is always understood 
to be a subset of H: 

DomJ := {rt G iL : J{u) < oo}. (2-2) 

Let fl be an open subset of W^. A function u G L^(n) whose gradient Du in the sense 
of distributions is a (vector-valued) Radon measure with finite total variation in El is 


u 


:={u,u)^/‘^, {u,v) := 
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called a function of bounded variation. The class of such functions will be denoted by 
The total variation of Du on Q is defined as 


|L)u|(0) = sup |y u divz dx : z G | 2 ;(x)| < 1 Vx G , (2.3) 

where for a vector v = {vi,... ,Vn) G we set |up := denoted by 

|Z)u|(n) or by \Du\. The map u —)■ |Du|(n) is L;^^^(17)-lower semicontinuous. BV{^) 
is a Banach space when endowed with the norm ||u|| := \u\ dx + |Z)tt|(n). 

Analogously one defines BViod,^) as the space of functions u G such that 

\Du\{U) < +00 for any [/ d 12. 


A measurable set i? C 12 is said to be of finite perimeter in 12 if (2.3) is finite when u is 
substituted with the characteristic function xe of E. The perimeter of i? in 12 is defined 
as P{E,Q) := \Dxe\{^)- We denote by and respectively, the n-dimensional 

Lebesgue measure and the (n — l)-dimensional Hausdorff measure in M”. 

In case ill is a set with finite perimeter, one can define its “reduced boundary” d*E as 
the set of points where there exists the limit 


DXE{Bp{x)) 

\DxE{Bp{x))\ 


ue{x) g5"-i 


and it has norm one. Then, one has the representation 

E>Xe = ve{x)\Dxe\ and \Dxe\ ^d*E. 


The first equality just follows from Besicovitch’s derivation theorem for Radon measures, 
while the second from a careful study of the reduced boundary, see Section 3 in [7]. In 
particular, it follows that for any open set A C 12, 

p{E,A) = r^-^{^*Er^A). 

If if C M" is a measurable set and x G ME, we define the (Lebesgue) upper density of 
LI at X by 

— _ , , |Llnii„(x)| 

D{E,x) := hmsup ' , , 

\^py^)\ 

while the density, when it exists, is simply the limit limp_ 5 .o \E n Rp(x)|/|i?p(x)|. Then 
it can also be shown that up to ^-negligible sets, the reduced boundary d*E of a set 
with finite perimeter coincides with the set of points where E has density exactly 1/2, 
and with the set of points where E has density neither 0 nor 1. 

Given u G BV{Q), we define 

u'''(x) := inf{2 ; D{{u > t},x) = 0} and u~{x) := sup{2 ; D{{u < t},x) = 0}. 

Then, we say that u is approximately continuous at x G 12 if and only if u'^{x) = u~{x). 
The set of points where u is not approximately continuous is called the singular set of 
u and denoted by Su- If "U is the characteristic xe of a set with finite perimeter, then 
d*E C and \ d*E) = 0. On the other hand, almost any level {u > 2} of a 

BV function has finite perimeter, and there holds the co-area formula 

/ + 00 

P{{u>t},n)dt. (2.4) 
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For a comprehensive treatment of functions of bounded variation, we refer the reader 
to [7]. 

2.2. Basic theory for J and Ix. In this section we assume that is a bounded, open 
subset of and that f,g satisfy (1.4), (1.5). We collect here some basic results about 
J and Ix, see e.g. [15, 4, 14]. We first characterize the domain of the functional J: the 
following useful fact is a special case of the general results established in [6]. 

Lemma 2.1. If u ^ then J{u) < co if and only if 

u G and f is \Du\-integrable, 

where \Du\ denotes the total variation measure assoeiated to the gradient measure Du. 
Moreover, J{u) equals the integral of f with respect to \Du\. 

We now recall a version of the coarea formula. The statement we give (in particular 
our convention for defining the sub/super levelsets Et) is more or less forced on us by the 
fact that if w G H and g 0 L^(fl), then X{w>t} 0 77 for t < 0. 

Lemma 2.2 (Coarea Formula). Assume that w G DomJ, and for t G M define 

Et = Et{w) := {x G n : sign(t)(r(;(x) — t) > 0}. (2.5) 

Then Et G DomJ for a.e. t, and 


J{w) = f f\D 
Jn 


w\ = 


'd*Et 


fdn 


n—l 


dt = 


J{XEt)dt. 


Proof. Since \w\ > jtj in Et, it is clear that Xst — and hence that XEt £ H for 

every nonzero t. Next, Lemma 2.1 implies that w G BVioc, and then the standard coarea 
formula implies that the set {x G ; w{x) > t} has locally finite perimieter for a.e. t, 
and that 


/ f\Dw\ = r 

J Q J 


fdn 


n—l 


dt. 


jQnd^{w>t} J 

Then the thesis follows from the facts that d*{w > t} DU = d*{w < t} CiU for every t, 
and that d*{w <t}nfl = d*{w > t} for every t such that the level set {w = 0} has 
measure zero, which holds for a.e. t. □ 


Similarly, we will sometimes use a “layer-cake decomposition” in the form 


(p{x)q{w{x),x) dx = 


sign(t) / ip{x)qt{t,x) dx dt 
JEt 


( 2 . 6 ) 


where w G H and ip G Li/g, and q{t,x) is a Lipschitz function (often q{t,x) = t) such 
that q{0, x) = 0. Indeed, in this case 

[ <p{x)q{w{x), x) dx = 


/ / X{o<t<w{x)}qt{t,x)dtdx- 

J {ui>0} J 

and then (2.6) follows by Fubini’s Theorem. 




/ X{o>t>w{x)}qt{t,x)dtdx, 
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Notice now that for a general convex positively 1-homogeneous function J on a Hilbert 
space H we have, for any u G Dom J, 

7] G dJ{u) J{v) — J{u) > {r],v — u) for any v G Dom J, (2-7) 

where (•, • ) denotes the inner product in H. Testing with v = Xu for A > 0, we obtain 

T]£dJ{u) J{u) = {r],u). (2.8) 

We deduce from (2.7) the equivalent characterization 

rj G dJ{u) J{u) = {rj,u) and (r?, u) < J{v) for any v G DomJ. (2.9) 

In particular, the convex conjugate function J* corresponds to the indicatrix function 
J*{r]) = 0 if snp{r],v) < J{v), and J*{ri) = +oo otherwise. 

The following lemma, a small modification of results of Anzellotti [8], will be needed 
for a characterization of the subgradient dJ. 

Lemma 2.3. For any v G DomJ C and ^ G L°°(H;M"') such that div(/^) G Li/g, we 
define a distribution (/^ • Dv) by 

[{f^-Dv)(f:=-[ [ div{f()v(f, for(peC^{n) (2.10) 

Jfi Ju Jq 

Then (/^ • Dv) is well-defined and satisfies 


in 


iffi ■ Dv)ip 


< llv^lloo IlCIloo / f\Dv\ 
Ja 


if supp{(p) C A open, 


( 2 . 11 ) 


and hence is a Radon measure with finite total variation. 


Proof. This is proved in [8, Theorem 1.5] if / = = 1. In the present situation, our stand¬ 

ing assumptions (1.4), (1.5) and the fact that div{ffi) G guarantee that the integrals 
on the right-hand side of (2.10) are well-defined. Then (2.11) is clear if v is smooth, and 
the general case follows, as in [8], by a density argument. This is straightforward since 
we only require approximation, in the interior of H, of u G Dom J by smooth functions, 
and this is clearly possible in view of (1.4), (1.5). □ 


Lemma 2.4. // n G DomJ, then rj G dJ{u) if and only if there exists a vector field 
f, G L°°{Q;MP) such that ||^||oo < 1? 


/ ifC-Dv) = {r},v) 
Jn 



for all V G Dom J, 


( 2 . 12 ) 


and 

[ iff-Du) = J{u) = [ f\Du\ . (2.13) 

Jn Jn 

In particular, dw{ff) G L^g, so that (/^ • Dv) makes sense for all v G DomJ. 

Moreover, if (2.12), (2.13) hold and Du is an integrable vector field, then ^ o.e. 

in supp (Du) n H. 
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Notice that, by taking v G C DomJ, we get that 

(2.12) ^ div (/0 = -V 9 in 

Formally, (2.12) also implies the boundary condition 

/ vf^ ■ V = 0 for all p G Dom J C Lil. 

Jdn 

One cannot immediately conclude that f^-iy = 0on dQ, even formally or in a weak sense, 
without further assumptions on g, since membership in H may force every v G Dom J to 
vanish to high order on sone subset of Ski. On the other hand, the boundary condition 
/^ • V = 0 trivially holds on portions of dVt where / vanishes. 


Proof. Fix rj G dJ{u), and define the Banach space 


Ai = Aif := {M’^-valued measures /2 on 11 : ||/l||x^ := / /d|/f| < oo}. 

Jn 

If u G Dom J, then Dv G Aij, and 

\{t],v)\ < J{v) = \\Dv\\Mf ■ (2.14) 

Thus the linear functional Dv i— {g, v) extends, by the Hahn-Banach Theorem, to a 
linear functional /? : A4/ —)• M satisfying 


/3{Dv) = {g,v) 



for V G Dom J 


(2.15) 


and 

|/3(m)| < f f d\il\ for all /i G A4/. (2-16) 

Jn 

Let /So denote the restriction of (3 to Lj.(0;]R”). From (2.16) we deduce that there exists 
a vector field with ||^||oo < I, such that 


Poi'ip) = Pi'if) = j^f{x)f{x) ■ i>{x) for all V) G L}(0;M”). (2.17) 

By combining (2.15) and (2.17), we see that if u G (^“(ll) C DomJ, then the distri¬ 
bution div(/^) satisfies 

^ (^^div(/^)^ (^u) :=J^f^-Dv = /3o{Dv) < C\\v\\h = \\y/gv\\L 2 (^<^y 

Since this holds for all v G C')T(ll), and since g is bounded below on compact subsets of 
n, it follows that div(/,f) G As a result, l3{Dv) = /Qdiv(/^)u for all compactly 

supported v G Dom J. 
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Next, for v G DomJ and ip G note that D{vp) = vVp + pDv, from which it 

follows that vp G Dom J. Thus, by the linearity of /3, 


f3{pDv) = /3{D{vp)) - l3o{v Vp) 

= - [ div{f()vp- [ fi-Vpv 

Jfi Jfi 

= [ {fi-Dv)p. 

Jn 


Since (/^ • Dv) is a finite measure, and in view of (2.16), this identity continues to hold 
for p G C{Q). In particular, we may set p = 1 to find that 

/ (/^ • Dv) = (3{Dv) = (?7, v) for v G Dom J. 

Jn 

Thus (2.12) holds, and (2.13) follows from (2.8). 

Conversely, if (2.12) and (2.13) are satisfied, then it follows from (2.11) that 

{r],v) < J{v) for all v G DomJ, = J{u). 

As noted in (2.9), this implies that rj G dJ{u). 

Finally, if u G then • Du G L^(n). So in this case (2.12), (2.13) imply that 

[ f {\Du\ -Du-0= f f {\Du\ -Du-0 = 0 

JQ ^supp(Dii) 


Llpp(Dli) 

which, since ||.^||oo < implies the final conclusion of the lemma. 
A consequence of (2.13) is the following 


□ 


Lemma 2.5. Assume the hypotheses of Lemma 2.f, and define Et = Efiu) as in (2.5). 
Then rj G dJ{u) if and only if sigiii{t)rj G dJ{xEt) for a.e. t G M. 

Proof. Assume that p G dJ{u) and fix ^ as in Lemma 2.4. Then, denoting by XEt the 
characteristic function of Et, Lemmas 2.2 (the coarea formula) and 2.4 imply 


/ J{XEt)dt = J{u) = [ iff,-Du) 
Jm. j D, 


(2.18) 


On the other hand, fixing a sequence of compactly supported smooth functions Pk Xn, 
we have 

[ iff. ■ Du) = lim f (/^ • Du)pk 
Jn ^ Jn 

= lim [ sign(f) [ {ff ■ DxEfiPk dt = [ sign(t) [ {fODxEt)dt. 

^ Jr Jn Jr Jn 

The second equality follows from (2.10) and (2.6), and the last equality follows from 
(2.11), (2.18), and the Dominated Convergence Theorem. Combining these identities 
yields 

[ dt(j{xEt) -sign(t) [ iff-DxEt)'] =0. 
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Since |^| < 1, by (1.1) we have, for any t G M, J{xEt) > sign(f) /^(/^ • DxEt)-, and hence 
for (2.18) to hold we must have, for a.e. t G M, 

J{XEt) = [ (/(sign(t)O • DxEt) ■ (2.19) 

Jn 

In view of (2.12), (2.13), this states that sign(t)ry G dJ{xEt)- The proof of the converse 
is very similar. □ 


2.3. Properties of minimizers of Ix. We have the following 

Proposition 2.6. Given f,g satisfying (1.4), (1.5), A > 0 and uq G H, there exists a 
unique minimizer w G Dom J of Ix. This minimizer is characterized by 

X{uo — w) G dJ{w). (2.20) 

In addition, if uq G L°°(n), then 

ess inf Uq < ess inf w < ess sup w < ess sup uq . (2.21) 

n Q n n 

Proof. Note that /a(0) = |||no|P < oo, so Domix is nonempty. The functional Ix is a 
coercive, convex, weakly lower semicontinuous functional on H, so that the existence of 
a minimizer w G Dom J is obtained by direct methods in the calculus of variations (see 
[7]). Uniqueness is a consequence of the strict convexity of the functional Ix- 

If we temporarily write Q{u) = ^||n — uqIP) then the Euler-Lagrange equation for the 
minimizer of Ix is 

0 G dlx{w) = d{J + XQ){w) = dJ{w) + XdQ{w), 

where the final equality follows from standard convex analysis considerations, see for ex¬ 
ample [17], Proposition 1.5.6. Upon rewriting this, we obtain the necessary and sufficient 
condition (2.20). 

The uniform bound (2.21) follows by noticing that, for any v G 
/a( ess infuoVuAess snpuo) < Ix{v). 


□ 


2.4. Dual formulation. Equation (2.20) is equivalent to 

w G dJ* {X{uq — w )), 

where J*, the convex conjugate to J, is given by 

J*(y) — / ^ ^ 

^ 1 +00 otherwise, 


( 2 . 22 ) 

(2.23) 


and 

||u||* = sup {{v,u) : n G H, J{u) < 1} . (2.24) 

We may view the norm || • ||* related to the dual convex function J* as a weighted 
norm. Very much as in the poof of Lemma 2.4, the Hahn-Banach theorem implies the 
variational characterization 
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The following variant of a standard fact illustrates the impact of the (non)integrability 
of g. 

Lemma 2.7. Assume that uq G H = Lg, let w he the minimizer of I\, and define 


a — dst: 


f g(x)uo(x)dx/J g(x}dx ifgGL^{^l) 
0 ifg^L^{n). 


( 2 . 26 ) 


Then 

w is constant <^=i> w = <^=i> A < Ac = -n-n- . (2.27) 

||uo — a* II* 

Moreover, if g G L^{Q,), then 

j g{uo -w) = 0. (2.28) 

The minimizer w does not need to satisfy (2.28) ii g ^ (when the integral makes 
sense, which is not guaranteed). This can be seen from the concrete problems considered 
in Section 5, in which uq is positive and belongs to but ||tto||* < oo, so that w = 0 is 
the minimizer for sufficiently small values of A. 

If n > 3, there exist uq G H such that ||wo — a*||* = +oo, and hence w is nonconstant 
for all A > 0. In the unweighted case, this is essentially equivalent to the statement that 
does not embed in for n > 3, and the presence of weights satisfying (1.4), (1.5) 

does not change things. Whether or not such uo exist in 2-dimensional domains depends 
on the behavior of / and g near d^l. 


Proof. It follows from (2.22) and (2.23) that 

A||uo — Ml* ^ 1) (2.29) 

since dJ*{X{uo — w)) = 0 if J*(A(uo ~ 'w)) = +oo. 

If g is integrable, then every constant function c belongs to Lg. Then (2.29) and the 
definition (2.24) of the || • ||* norm imply that (uq — w,c) < ||tt||* for every c, and hence 
that {uq — w,1) = 0, which is (2.28). It then follows from (2.28) that if w is constant, 
then the constant must equal a*. 

On the other hand, if g ^ L^, then no nonzero constant belongs to H, and so the only 
possible constant minimizer is a* = 0. 

In either case, the definition of Ac and the characterization (2.23) of J* imply that if 
A > Ac then J*(A(tto — a*)) = -|-oo. Then d.J*{X{uo — a*)) = 0, and w = does not 
satisfy (2.22). 

On the other hand, if A < Ac then J*{X{uo — a*)) = 0, and it follows that a* G 
dJ*{X{uQ — afij). This is clear if a* = 0, since J* attains its global minimum at A(tto —a*). 
If a* 0, then necessarily g G L^(n), and then one can check that every constant function 
belongs to dJ*{X{u — a*)). Indeed, for any v G H and any a G M 

J*{X{uo — w)) + {a,v — X{uo — a*)) = (a, v) < J*(u). 


□ 
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Remark 2.8. Using the language of differential forms, and fixing a 1-form oiq such that 
duo = 0 (i.e. loq is closed), (/wq) at = 0 on dO, \ K and v = d*{f ■ wq), we may reformulate 
the variational problem (2.25) as follows 

, G C“(fl;A2(M^)), V’A = 0 on . (2.30) 

OO J 

If one seeks a 2-form il) such that {d*ip/f) -|-a;o minimizes the L°° norm with respect to all 
compactly supported perturbations of i/’; then the resulting problem may be interpreted 
as a sort of (weighted) vector-valued analog of the Aronsson problem. It would be inter¬ 
esting to give further equivalent descriptions of the problem (2.30) in a game-theoretical 
framework, as the tug-of-war description of Aronsson problem (see e.g. [11], [24]). 


= inf 


d*^p 


f 


+ Wo 


2.5. Level set formulation. We shall say that a set i? C is admissible if xe £ Dom J. 
For an admissible FI C 11, we set 

P{E) = J{xe) = [ f{x)d'}r~^, A{E, t) = sign{t) [ {uq - t)g{x) dx . 

Ja*Enn Je 

and 

h,t{E) = P{E) - XA{E,t) 

Note that A{E,t) = sign(t)(uo — tXE,XE), which makes sense for an admissible set E. 
Lemma 2.9. If w £ DomJ, there holds 

h{w)-h{0)= [ Ix,t{Et)dt. 

Jr 


Proof. The thesis follows from the coarea formula (Lemma 2.2) and the “layer-cake de¬ 
composition”, see (2.6), applied to 

^ |(w - uof - j |uo = g{x)q{w{x), x) dx for q(t, x) = - tuo{x), 

where Et = Et{w) is defined in (2.5). □ 


Proposition 2.10. A function w £ DomJ minimizes I\ if and only if, for a.e. t, the 
set Et = Et{w) minimizes Ix^t among all admissible E C II. 

Proof. It is immediate from Lemma 2.9 that if Et minimizes Ix^t for a.e. t, then w 
minimizes Ix. 

Conversely, notice first that from the Euler-Lagrange equation (2.20) for Ix and Lemma 
2.5 we have sign(f)A(uo — w) £ dJ{xEt) fo^' a-e- 1 G M. Hence, if xe £ DomJ we have 

sign(t) [ X{xF - XEt){uo - w)gdx < J{xe) - J{XEt) = Pi^) - P{Et)- 
Jn 
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Moreover, (2.5) implies that sign(t)(t — w){xf — XEt) ^ 0; so we deduce that 


XA{F, t) - XA{Et,t) = sign{t) f X{xf - XEt){uf) - t)g dx 

Jn 

< sign(t) / X{xE - XEt){uo - w)gdx < P{F) - P{Et), (2.31) 

Jn 

and the conclusion follows. □ 

Remark 2.11. From the proceeding proof we deduce that if F minimizes /a,* then 
necessarily we must have equality in (2.31) and hence the minimizer w satishes 

{t — w{x)) ■ {xf{x) — XEti^)) = 0 for a-.e. x G Q,. 

This implies that Ft C F C := {x : sign{t){w{x) — t) > 0}. The sets Ft and E^ are 
called respectively the minimal and the maximal solution of the minimization problem 
min/A,i- In particular, since for a.e. t G M the level set {tc(x) = t} has zero measure. Ft 
is the unique minimizer of Ix,t for a.e. t G M. 

Remark 2.12. By (2.20) and (2.9) we deduce that 

{X{uo - w),xf) = [ X{uo-w)g{x)dx < [ f{x)dn"'~^ = P{F) (2.32) 

Jf Jd*F 

whenever xf G DomJ, with equality for F = Ft, for a.e. t G M, i.e. 

sign(t) f X{uo — w)g{x) dx = j f{x)dT-F~^ = P{Et). (2.33) 

JEt Jd*Et 

One may wonder if for F satisfying (2.33) we necessarily have that F minimizes /a,!- 
Since sign(t)(t — w)xEt < 0 it follows from (2.33) that for a.e. t G M, 


min/A,i = h,tiEt) = P{Et) - XA{Et,t) = sign(t) [ X{t - w)g{x) dx < 0 . 

JEt 


Moreover, min/A,t = 0 if and only if t > ess supj^ rc > 0 or t < ess inf^ tc < 0. 

Remark 2.13. Since I\^t{E) < Ix^s{E) for every E, when 0 < t < s, we have 

min/A^t < min/A^s for any 0 < t < s, 

Moreover, still assuming s > t > 0, we have Eg C Ft and 

t{s) — i{t) ^ = ^ — f it — w)g{x) dx — f {t — s)g{x) dx , 

'J Et\Es Es 


where z(r) = min/A.r) so that 


lim 

s—>-t+ 


s — t 



g{x)dx, 


lim 

t—¥S 


^js) - i{t) 

s — t 



g{x)dx, 


(2.34) 


(2.35) 


where FI® = n^<s= {rc > s}. Hence the map i possesses bounded left and right 
derivatives at every point. Flat level sets of w correspond to points of non-differentiability 
for i. Similar considerations apply to t < s < 0, with signs reversed. 
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3. Extremal level sets of the minimizer 

When (2.27) is violated, the minimizer w is nontrivial, and hence one may wonder 
whether w has “flat spots” - that is, whether any level sets of w have positive measure, 
as happens in the classical obstacle problem when the constraint is saturated or also for 
related models (see e.g. [13]). 

The answer to this question requires more assumptions on the boundary behaviour of 
/ and g with respect to those that we have made up to now. In particular, we will still 
assume (1.4) and (1.5), and we will also require that 

fix) > Cp^ix), |no(x)| gix) < Cp(x)-^, (3.1) 

for some positive integer m G N and some 6 G (0,1] to be specified below, where 

p(x) := dist(x, (912) for x G 12. (3-2) 

We continue to allow non-integrable ^’s, but we require that the possible blow-up of g 
near <912 is compensated by the decay of uq. 

Proposition 3.1. Assume that Q. <Z MA is a bounded open conneeted set with smooth 
boundary, and let w minimize I\. 

If f,g,uo satisfy (3.1) for some m G N and 0 G (^^ 7 ^q(^,l], then there exists some 
constant a > 0, depending on /, g, uq, n, 12, A, such that 

CI^iEt) > a for every t sueh that C^iEt) > 0. (3-3) 

In partieular, w G L°°(12), and if either of t^ := ess sup^w^ is positive, then the 
corresponding level set E^ := {x G 12 : wix) = ±2^} has measure at least a. 

Finally, if (3.1) holds with 9 = , then w G L°°(12). 

Remark 3.2. Proposition 3.1 remains valid if there exists h : —)• M such that 

12 = {(xi,x') G M X : 0 < xi < h{x)}, |V/i(x')| > 0 where /i(x') = 0, 

and in addition / > 0 on {x G c212 : xi =0}. Indeed, in the proof of the proposition, the 
smoothness of (912 is only used in invoking the weighted relative isoperimetric inequality 
(see Proposition 3.4 below), and this result is still valid under the hypotheses described 
here. 

Remark 3.3. If uq > 0, so that in > 0 and t~ = 0, then Proposition 3.1 does not in 
general say whether the level set {w = ess inf^rc} is flat. 

However, if gix) < (which in particular implies that g is integrable) and uq G 

L°° , then both extremal level sets of w are always flat. Indeed, the hypotheses of the 
proposition are still satisfied by 

Ixiu) := Jiu) + ^ [ g\u-uo\'^, ioi uq := uq - c 

2 Jn 

for any constant c. Moreover, if w minimizes Ix, then w := w — c clearly minimizes lx- By 
a suitable choice of c, one can arrange that fio < 0 in 12, and hence that ess inf^ w < 0. 
It follows that the level set {w = ess infn w} = {w = ess inf^ re} is flat. 
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On the other hand, in the model case discussed in the next section, the results in [2] 
imply that {zc = ess infn zc} may have either positive or zero measure, depending on the 
choice of the parameters. Note that in this case, we have g 0 

The main tool in the proof of Proposition 3.1 is the following weighted relative isoperi- 
metric inequality. 


Proposition 3.4. Assume that Q, C is a bounded open conneeted set with smooth 
boundary. 

Then for every m gN, then there exists a constant C = C{Q,,m,n) such that for every 
set E C id of locally finite perimeter, 


mm 


( [ p^{x)dx, 

f p^{x)dx\ <C{Q,m,n) 

f p^{x)dn^-^ 

\Je 

Jn\E J 

Jd*Er\n 


n-\-m 
n + m —1 


(3.4) 


The proof is deferred to the Appendix. For now, we accept the result and use it to 
prove that positive and negative extremal level sets of the minimizer are flat. 


Proof of Proposition 3.1. We first claim that there exists a constant ai > 0 such that 

if A{E,t) < oi and 0 < f dx < - f p'^ dx, then > 0. (3-5) 

Je ‘d Jii 

Note that under these conditions, our hypotheses on / and the relative isoperimetric 
inequality (3.4) imply that 

n-\-m 

f if 1 \ '^+^“1 n+m 

0< / p^dx <c{ fdn^-^) =CP{E)^+^k (3.6) 

Je \Jd*Enn ) 

For simplicity, we will prove (3.5) for t > 0; the argument for t < 0 is essentially 
the same. First, we may assume that A(E,t) > 0, since if A{E,t) < 0, then Ix^fiE) > 
P{E) > 0 by (3.6). 

Next, for an admissible E with A(E, t) > 0, we compute 


A{E,f)= / {uQ — f)gdx< / {uq — f)'^g dx < C / p' 
Je Je Je 


Pm ( up - ty 
uo 

using our assumption (3.1) on g. Thus by Holder’s inequality and (3.6), 


dx 


for 


A{E,t) <c(^J p™ dx 
6(n + m) 

0- = —-r > 1 , 

n + m — 1 


It follows that 


hAE) > 


' (up - 1) 
In V ^^0 


c{t) = C 

A{E,tA^/'^ 


— f\+\ i-e 


1-0 


dx < c{t) P{Ey 


(zzo - tA\ 


1-0 


c{t) 


n\ Up 


- XA{E,t). 


dx 


Since c{t) < CC"'{QA ^ for all t, is is easy to find oi > 0 such that Ix,t{E) > 0 if 
0 < A{E,t) < ai, proving (3.5). 
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Next, fix a > 0 such that if C^{E) < a, then the hypotheses of (3.5) are satisfied. If 
0 < t we deduce from (2.34) that < 0 and then from (3.5) that C^{Et) > a, 

proving (3.3). It easily follows that w G L°° and that C^{E'^) = £”'(no<t<t+£^t) > a. 
Finally, ii 0 = , then the above arguments show that if A(Et,t) > 0, then 


By the Dominated Convergence Theorem, there exists some T > 0 such that c{t) < ^ 
whenever t >T. As a result, if t > T, then A{Et, t) < 0, which implies that I\^t{Et) > 0, 
and hence (again by Remark 2.12) that t > = ess supj^tc. This proves that < T. 

Since similar considerations apply to t~ = ess sup^ w~, it follows that w G L°°{Q,). □ 


0 > It,\{Et) > A{Et,t) { 


4. BV AND BOUNDARY BEHAVIOR FOR SMOOTH WEIGHTS 

In this section we impose further decay and smoothness conditions on the weight /, as 
well as compatible decay conditions on g, involving the minimizer w as well as the data 
Uq. 


Lemma 4.1. Assume that f = for some positive integer m and oj G (7^(11) such that 

iVtcp + tc > Cl > 0 in Q. (4.1) 

Let w minimize I\ and assume that there exists some eonstant Ci such that 

\g{w — uo)\ < Cif a.e. in LI. (4.2) 

Then u G BV{Ll), and there exists some constant C, depending on k, ||w||c 2 )Ci and Ci, 
such that 


\Dw\ < CJ{w) 


and LL 


n—1 


{d*Et nn)<C f fdn^-^ for a.e. 


Assumption (4.2) is in general hard to verify, since it requires knowing something about 
the minimizer w. However, 

• it is satisfied A g < Cf and uq G T°°. 

• If (/ does not decay at dLl, and in particular if 5 0 L^(D), it can also be proved 
to hold for certain uq. Indeed, we prove in Section 5 that this is the case for the 
problems related to Bose-Einstein condensation. 

The proof is an adaptation to our setting of an argument from Section 3 of [23], the 
main difference being simply that we use the first variation formula and an approximation 
argument in place of the direct ODE arguments of [23]. 

Proof. Consider V G and let 

4>^(a:) := X + hV{x) 

Then there exists /iq > 0 such that if \h\ < ho, then is a diffeomorphism of Ll onto 
itself which fixes a neighborhood of dLl. 

Recall from Proposition 2.10 that Et = Et{w) minimizes = P {')~among 

all admissible sets E, that is, sets such that xf £ Dom J. Consider any t for which this 
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holds, and for \h\ < ho let := Since is C^, standard theory implies that Ej^ 

is a set of locally finite perimeter in f]. Moreover, E^ coincides with Et in a neighborhood 
of (9n, and thus is easily seen to be admissible. Then the minimality of E\^t implies that 


d 

~ ~dh 


hAE^ 


h=0 


ld*Et 


f{dij-EiUj)VA+VAxi\dn'^ -Asign(t)/ g{uo - t)V ■ v dU 

J Jd*Et 


n—1 


where v{x) denotes the approximate outer unit normal vector to Et Sil x ^ d*Et. Indeed, 
the first integral on the right-hand side is the first variation of the weighted perimeter 
P(-), see for example [21], chapter 17, and the second integral on the right-hand side is 
the first variation of 

Let now V = q{uj)Vuj, for some smooth function q such that 

0 < g < 1, q > 0, (7 = 0 in a neighborhood of 0. 

Since q o u, and hence V, are compactly supported, we deduce that 

f q{Ld)muj"^~^\Vuj\‘^d7i'^~^ = — f — p (g) p) ; (Vw (g) Vcu) 

Jd*Et Jd*Et 

- f u^q{f){I-u®u) 

Jd*Et 

-Asign(t)/' q{f)g{uo-t)Vuj ■uddE~^. 

Jd*Et 

For a.e. t, it follows from (4.2) that \g{uo — t)\ < C\f = Cioj^ almost everywhere on 
d*Et. Using this and the fact that q' > 0, we deduce that 

jn—1 


'dn 


j q{(jj)muj'^ [ u}"^q{uj){I — u 0 u) ■dV. 

Jd*Et JdtEt 

- A sign(f) / q{uj)g{uo -t)v-Vf 
Jd*Et 

< ||u;||c2(1 + AC'i) / 

Jd,Et 

Letting q A Xir+ > we conclude that 

f oj^-^\Vio\An^-^ < c f oj'^d'H^-K 

Jd*Etnn Jd*Etnn 

for C = C'(||a;||f 72 , Cl, A), for a.e. t. Since jVwp -|-w > ci > 0 in U, it follows that 


^yn-l < 1 


I d* Etna. 


[ + ^)d'H'^~^ <C [ 

Jd*Etnn Jd 


U! 


d*Etnn 


Cl Jd*Etr\n 

Repeating essentially same argument, with V replaced successively by q{uj)uj~^Vuj for 
j = 1,..., m — 1, we find that 

f dn^-^ <C [ dTL^-^, 0<j<m-l 

Jd*Etnn Jd*Etnn 
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for a.e. t. Putting together these estimates, we conclude that < CP{Et) for 

a.e. t. The estimate \Dw\ < CJ{w) then follows from the coarea formula. □ 

The next result is due to Montero and Stephens (see Section 4 of [23]) in a slightly 
different setting. For this we specialize to n = 2, and we impose still stronger assumptions 
on f,g,uo. 

Lemma 4.2. Assume that n = 2, that the hypotheses of Lemma 4-1 are satisfied, and 
that in addition 

g G uo G TT1’°°(S4). 

Then for a.e. t, dEt is a union of curves meeting orthogonally dLl . 

We recall the argument from [23]; this will be needed later when we argue that the 
result still applies to certain models of Bose-Einstein condensates for which g 0 L^{Q,). 

Proof. Fix some t such that Et minimizes and Pfi{d*Et) < oo and such that for 
(3 := g{uo — t), we have |/3| < Cif almost everywhere on d*Et. (These conditions 
are satisfied a.e. by Proposition 2.10, Lemma 4.1 and (4.2).) The minimality of Et, 
together with standard regularity results, implies that dEt n 11 is a countable union of 
nonintersecting C^-curves, each of which must have finite arclength. Let 7 : (a, 6) —)• 11 
be an arclength parametrization of one such curve. The Euler-Lagrange equation for Ix^t 
may be written (along 7 ) in the form 

- (/( 7 ('S)) 7 '('S))' + V/( 7 (s)) - A/3(7 (s))p(s) = 0 

where v(s) is the outer unit normal to dEt at 7 ( 5 ). Suppose that 7 ( 5 ) approaches <911 as 
s a or s b; we consider the former case for concreteness, and we assume that a = 0 . 
Then by integrating the above equation from 0 to s, taking the inner product with p(s), 
and rewriting, we find (since / = 0 on 911) that 

0 = -/(7('S))7'('S) • p(s) 
s 

= ^{s) ■ V/(7(s)) + - / ■ \ (V/( 7 (l)) - V/( 7 (s))) - A/ 3 ( 7 (l))v(l)j dt . 

s Jo 

Using this and hypotheses (4.1), (4.2) one can verify that u{s) • Vuj(j(s)) —)■ 0 as s \ a, 
which implies that the curve meets 911 orthogonally. 

□ 

5. Applications to quantum fluids 

As mentioned in the Introduction, in certain limits the velocity field in rotationally 
forced Bose-Einstein condensate is described by a particular instance of the functional Lx 
from ( 1 . 2 ), where A = 2 and /, g, H have the form 

H := {(r,z) G (0,oo) X M : pj,^(r, 2 :) > 0}, f = Ptf^ g{r, z) = (5.1) 

Here Pj,p ; —)• M is a smooth function known as the “Thomas-Fermi density” and 

determined, in physical experiments, by the form of a potential that is used to confine 
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the condensate. We will assume that 

Ptf + > 0 on 90, is 0 is bounded. (5.2) 

The rotational forcing is encoded in the function uq, which has the form 

UQ{r,z) = ar‘^ (5.3) 

where a > 0 is a parameter. 

The following results about this functional are mostly immediate consequences of the 
general theory developed above. 


Theorem 5.1. Let /"(•) denote Ix{ ) for A = 2 and f,g,uo as in (5.1), (5.3), where 
satisfies (5.2) and either 

^Ptf z) < 0 in Q (5.4) 

or 
or 


Ll is a compact subset of (0, oo) x M. 


(5.5) 


Then for every a > 0, there exists a unique minimizer Wa G DomJ of 

(1) Wa is the unique element of H = L'^ satisfying the equivalent conditions: 

• the Euler-Lagrange inclusion (2.20); 

• the dual formulation (2.22); 

• the level set formulation from Proposition 2.10. 

(2) Wa is constant (corresponding to a vortex-free condensate) if and only if 


0 if (5.4) holds 

J^gr'^drdz/ J^gdrdz if (5.5) holds 

where the || • ||* norm is defined in (2.24). 

(3) If we write Ifi{E) := f dPL^ ~ ^ Ie ^(“^^ “ 

||R^a||L°° = inf{t > 0 : If{E) > 0 for all admissible E C O}. (5.6) 


a < 


Ir^ — 


for c* = 


(4) The level set {x G LI : Wa = maxniCo,} has positive measure. 

(5) If (5.5) holds, then the level set {x G LI : Wa = inf^tCc} has positive measure. 

(6) If (5.4) holds, then 0 < Wa < uq. As a result, Wa = = 0 on 90 n {r = 0}. 

(7) Wa G BV{Ll), and = 0 on 90 n {r > 0} in the sense that a.e. level set of 

Wa meets dLl orthogonally. 

(8) If (5.4) holds, and if E) denotes a connected component of Et for some t > 0, 

then El is simply connected and dE] n 90 0. 


Remark 5.2. Both the hypotheses (5.4), (5.5) are motivated by physical experiments. 

Indeed, assumption (5.4) covers any convex trapping potential that is symmetric about 

the vertical axis. In particuar, this includes the model case of a quadratic trapping 

2 2 

potential, corresponidng to Py^(r, z) = po ~ ^ ~ for positive constants po,a, 6. Note 
that when (5.4) holds, necessarily g 0 L^(0). 

The second case (5.5) is relevant to toroidal Bose-Einstein condensates. 
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Remark 5.3. Recall that the level set formulation states that for a.e. t > 0, the level set 
Et = {re > t} of Wa minimizes the functional /“ appearing in conclusion (3). It follows 
that dEf is analytic in 17 and satisfies the equation 

div{pTpU) = PppK + Vppp ■ V = ‘Ippp (^ar - (5.7) 

where p and k denote respectively the exterior unit normal and the curvature of dEf. 
(This equation appears already in a different form in the proof of Lemma 4.2 above.) 

Proof of Theorem 5.1. In case (5.5), / and g are smooth in Cl, and both are comparable 
to the distance p{-) to dCl, so all the hypotheses of every result in Sections 2-4 are 
satisfied (with m = 1 for Section 3), and all the conclusions above follow, after noting 
that conclusion (3) is a direct consequence of the level set formulation. 

If (5.4) holds, then necessarily p^p > 0 on {(r, z) £ dCl : r = 0}, and the form (5.1) of g 
implies that g 0 L^(I7). In this case, the results of Section 2 are still available, and these 
imply conclusions (1) - (3) of the theorem. Moreover, it is clear that \uo\g = arp^p < 
Cppp = Cf, so the hypotheses of Proposition 3.1 (see also Remark 3.2) are satisfied, and 
it follows that conclusion (4) holds. 

We next prove conclusion (6). The inequality Wa > 0 follows by noticing that /“(O V 
w) < I^{w) for all w £ DomJ. The inequality Wa < is equivalent to the inclusions 

Et C Hf for a.e. f > 0, (5.8) 

where we set = {ar^ > t}. To prove (5.8) it is sufficient to show that 

It{E n Hf) < It{E) for admissible any E C Cl. (5.9) 

In fact, recalling that < 0 and < t in E\ Ht, we have 

[ Ptf — [ Ptf f —f)< f (ar^ — t) 

JdpEnH^) Jd*E Je\H^ f Je ^ 

for any admissible E, which gives (5.9). 

It follows that the hypotheses of Lemma 4.1 are satisfied (in particular condition (4.2) 
is a consequence of conclusion (6)), and hence that w £ BV{Cl). 

Next we prove that = 0 on dCl n {r > 0}. It suffices to show that for a.e. f > 0, 
the conclusions of Lemma 4.2 hold for d*Et. The lemma as stated does not immediately 
apply here, since it assumes that g £ IT^’°°(I7). But the proof of Lemma 4.2 only requires 
g to be Lipschitz near d*Et, and this is the case here by (5.8), so the conclusion follows. 

Finally, if gy < 0 in 17, then we have dEt H dCl 0, since otherwise one could move 
Et in the r-direction and decrease its energy. Moreover, since — 7 < 0 in Fit, it also 
follows that Et “cannot have holes”, that is, each connected component of Et is simply 
connected. □ 

We now illustrate a consequence of conclusion (5.6) above. 

Lemma 5.4. Assume the hypotheses of Theorem 5.1 and assume that ^ < 0 in Cl. For 
t > 0, if there exists a 1-form Cf such that 

(1) dC^ = —2ppp{ar — ^)dr A dz in Cl 

(2) C* = 0 on on {{r, z) £ dCl : r > 0} 
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(3) — |C*| > 0 in = {(r, z) £ VL : ar"^ > t} 

then ||wq,||oo < t. 

Proof. Properties of C* and Stokes Theorem imply that for any nonempty admissible E, 

if{E)= f pdn^+ f C> j {p-\cf\)dn^ >0, 

Jd*E Jd*E Jd*E 

where d*Et is oriented in the standard way in the integral The final inequality 

above uses conclusion ( 6 ) of Theorem 5.1. It follows that Et, as the minimizer of 
must be the empty set, and hence that ||tCQ,||oo < i by (5.6). □ 


On the other hand, (5.6) also implies lower bounds for ||tCQ,||Loo, such as 
if if ( 7 ) := n {r > 7 }^ < 0 for some 7 > 0 , then ||rcQ,||oo > i- 


(5.10) 


Example 5.5. We now consider the model case arising from a quadratic trapping po¬ 
tential: 

Ptf — Po - 2 ~ Po,CL,b positive constants. (5-11) 

In this case, a 1-form f satisfying the conditions of Lemma 5.4 is given by C* = 7 *(r, z) dz, 
where 


a a 


aR\ 


r j 


T] := — -t[2R In ( — ) - ) , where R = R{z) -.= [ pQ - — 


v2\ 2 


Since p* = 0 on {(r, z) G dPl : r > 0}, and ^ < 0 in we see that p* > 0 in H\ 


(5.12) 


Thus 




t(a) := inf |t G [ 0 , aa^] (p^,^ — p*)(r, z) > 0 in iLf |, (5.13) 

for p* as in (5.12). 

On the other hand, a direct computation shows that the function if (•) appearing in 
(5.10) is 




8b 


7,,ptanh-‘(^J-7Z“-7«,Z 


Z\ 86 


86 


where Z = Z{'y) := (po — Thus (5.10) becomes 

ll'R'olloo > 1 ( 0 :) := sup < i G [0, aa^] min if(7)<0> 

1 7e(0,a^/p(f) ) 


(5.14) 


for i^ displayed above. 

Both t{a) and t{a) are easy computed numerically with considerable accuracy. In the 
next section we report the results of some such computations, and we compare the results 
to numerical solutions of the minimization problem. These results show that the upper 
and lower bounds for UrCaHoo given by (5.13) and (5.14) are often surprisingly close to 
each other. 
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Remark 5.6. The t = 0 functional Iq{-) determines whether or not ||rCa||oo > 0, and 
also determines the curve around which vortices first appear, at the critical value for a 
at which Wa becomes nontrivial. Following Lemma 5.4 and (5.12), we have 

Iq{E) = j pdli^ + l-aa^ [ dz. 

Jd*E ^ Ja*E 

This functional is considered in [2], where it is shown that if 6 > a in (5.11), then for 
every a > 0, the minimizer is either 11 or the empty set, whereas if a > cq • 6, where cq is 
an explicit constant, and a is not too large, then 11 does not minimize Iq. 

6. Numerical results 

Since the functionals studied in this paper are weighted versions of the Rudin-Osher- 
Fatemi denoising model [25], we can compute minimizers by adapting techniques devel¬ 
oped in the image processing literature. In the computations reported below, we have 
followed a modified version of a primal-dual algorithm [26], delineated in Algorithm 1 
below, to compute minimizers of the functional from Theorem 5.1, in the model case 
of a quadratic trapping potential, in which Pj,p is given by (5.11). We recall that the 
algorithm numerically solves the Euler-Largrange equations (2.20), see also Lemma 2.4, 
which in the present case may be written formally as as 

-rdiv(^/9j,^|^)-F2pj,^(u-ar2) = 0. (6.1) 

These are augmented with boundary conditions from conclusions (6) and (7) of Theorem 
5.1. Since the problem is formally overdetermined on cAl n {r = 0}, we only enforce the 
Dirichlet condition u = 0 there. 


Algorithm 1: Modified Zhu-Chan primal-dual algorithm. 

Data: Given a and TOL (tolerance for stopping) 

Initialize: k = 0,u^ = ar‘^,fP = |^[[o| ; 

while \\u^ — u^“^lloo > TOL do 

Choose time discretization parameters ri, and T 2 ] 

(we follow [26] and choose ri = 0.2 -|- O.OSfc, T 2 = . 

pk+i ^ pk _ 

k+i ^ 

” max{ , 1} ’ 

y^k _ yk-l _|_ ^^^2 _ yk-^ _|_ ^ (pP^)^ ; 

Set u^{r = 0) = 0; 
k = k + 1; 

end 
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We performed computations for the lower and upper limits of ||rco||oo found in (5.13) 
and (5.14) above, and we compared these to ||tCo||oo for an approximation Wa, obtained 
using Algorithm 1, of the minimizer Wa- 

In Tables 1 and 2 below we show the results of these computations, for a range of values 
of the parameter a and for two different choices of the parameters a and b appearing in 
the definition (5.11) of Pj,^. Numerical minimizers Wa for several choices of parameters 
are displayed in Figures 1 and 2. The computations show clearly that as the rotational 
forcing increases, the size of the vortex-free zone (the “flat spot”) decreases but does not 
vanish. 


a 

1(q) 

Oi 

||oO 

po ) 

Oi 

5 

0.140 

0.146 

0.21 

10 

0.321 

0.320 

0.39 

20 

0.485 

0.490 

0.54 

40 

0.6185 

0.620 

0.66 

80 

0.722 

0.723 

0.76 

160 

0.799 

0.800 

0.82 

320 

0.855 

0.856 

0.87 

640 

0.8967 

0.896 

0.91 

1280 

0.926 

0.925 

0.93 

2560 

0.9477 

0.946 

0.95 


Table 1. Numerical approximations of ||iCa||oo, and for various 
values of a. The experiments are performed for pancake-shaped regions 
i.e. po = 1, a = 1, and b = 0.5. 


Appendix A. A weighted relative isoperimetric inequality 
In this Appendix we present the proof of the weighted relative isoperimetric inequality. 

Proof of Proposition 3.4- Step 1. We first note that it suffices to show that there exists 
some Co > 0 such that 


/9™(x) dx < C 


ld*Er\n 


p^{x)dPL 


n—1 


n + m 
n+m —1 


if / /9™(x) dx < cq. (A.l) 
Je 


Indeed, for any 0 < cq < ^ /q P™", if we define 


ci(co, n) := inf I [ p^dn^-^ : CO < / P™ < J / / 

{Jd*Enn JE ^ Jn 
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a 

1(q) 

Oi 

W'^a ||cx) 

iM 

a 

20 

0.0350 

0.038 

0.05 

40 

0.0803 

0.085 

0.10 

80 

0.1211 

0.125 

0.14 

160 

0.1546 

0.158 

0.17 

320 

0.1802 

0.183 

0.19 

640 

0.1997 

0.201 

0.21 

1280 

0.2137 

0.214 

0.22 

2560 

0.2242 

0.224 

0.23 


Table 2. Numerical approximations of ||r(^a||oo) and for various 
values of a. The experiments are performed for cigar-shaped regions i.e. 
po = l,a = 0.5, and 6=1. 


then standard arguments imply that the infimum is attained. The fact that is smooth 
and connected implies that 

ci(co,0) > 0, 

and this fact and (A.l) together imply (3.4). 

Step 2. For 6 > 0 to be fixed below, we will write 

:= {x G ; p{x) < 6}, := {x G : p{x) > 6}, 

We will chose 6 small enough that is connected and is smooth, and moreover 
consists of a finite number of connected components, say 12* for i = 1,... ,k, one for each 
component of 512. (We will also impose a couple of additional smallness conditions on 6 
below, see (A.8), (A.9).) We will write 

E^ = EnQ^, Es = Enns, Ei = Enni, i = l,...,k. 

In view of (A.l), we may assume that 


p™ < ^ min 


P 


P 


/Di 


jQk 




P 


(A.2) 


Since 512'^ is smooth, it follows from (A.2) and the relative isoperimetric inequality in 12'^ 
(see for example [18] section 4.5.2) that 


Ie^ 


pm < ^ %'^-^[d*E^ n 12*') 


n/n—1 


< c 


VJd*Er\Q.^ 


p^dW 


n/n—1 


where all constants may depend on <5 as well as n, 12, m. If f^s p^ > k Ie P^^ then (A.l) 
follows easily, since then f^p"^ < C(f^s p™")^ for 6 := G (0,1). So we may 

assume that this inequality fails, i.e. that 



(A.3) 
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Once we assume this, in view of (A. 2 ) it is enough to prove the following weighted relative 
isoperimetric inequality in Oj for every i = 1 ,..., fe; 


1 


3™- 

J Ei ^ J J Ei 

Indeed, once this is know, then for E satisfying (A.2) and (A.3), we have 


if / < - / P^, then / p^<Q 

^-/o. .Iw.^ ^ -ld*Eir\Q.. 


p^dH^-^ 


n-\-m 


(A.4) 


P^<2^ P^<CY^ 


Er 


mjnjn-1 


'd*Eir\ni 


p^dn 


n+m 
n + m —1 


< 


^ E 


’d*Enni 


p^dTT-^ 


n+m— 1 


< C 


'a* EDO 




n+m 
n+m— 1 


Step 3. We now introduce auxiliary functions that will be used in the proof of (A.4). 
We will extend p to the complement of 0 by setting 


p{x) := 


dist(x,cI 0 ) if X G n 


—dist(x,cAI) if not. 

Thus the extended p is the signed distance to 50. We recall some standard facts about 
the signed distance function. First, because O is bounded and 50 smooth, p is smooth 
in 

Ns :={xeW^ : |/ 9 (x)| < 5}, 
for all sufficiently small (5 > 0 , and in this set, 

7 r(x) := X — p{x)Vp{x) is the (unique) closest point in 50 to x. 

In addition, 

Vp{x) = p( 7 r(x)), where p(y) = the inward unit normal to 50 at y G 50. 

In particular, 

|V/ 9 p = 1 , and by differentiating, V'^pVp = 0 . (A.5) 

For (x, y) G Ns X we will write 

p{x,y) := {p^{x) + = the (unsigned) distance from (x, y) to 50 x {0} 

O 5 := {(x,y) ; y(x,y) < 5}. 

Note also that for (x, y) G O 5 , the closest point in 50 x {0} to (x, y) is 7f(x, y) := (7r(x), 0). 
Now we define F : 0^ —)• O 5 by 

F(x,y) = 7 r(x)+ y(x,y)Vy(x) = x + [y(x, y) - y(x)]Vy(x). 

We note for future reference that F~^{x) is a m-dimensional sphere of radius p{x) for 
every x G O^; This can be seen by inspecting the definition of F, or it can be deduced 
from the formula 

F“^(x) = {( 7 r(x) + sV/j(x),y) : + |yp = p{xf], (A. 6 ) 

which is not hard to check. 
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Figure 1. minimizer Wa for a pancake-shaped condensate, with param¬ 
eters as in Table 1, for a = 20,40, 80 from top to bottom. 

Left: vertical axis is Waja. Note the change in vertical scale. 

Right; level sets of minimzers. 



























26 


P. ATHAVALE, R. L. JERRARD, M. NOVAGA, AND G. ORLANDI 



Figure 2. minimizer Wa for a cigar-shaped condensate, with parameters 
as in Table 2, for a = 20,40, 80 from top to bottom. 

Left: vertical axis is Waja. Note the change in vertical scale. 

Right; level sets of minimzers. 
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Step 4. It suffices to prove (A.4) for Ei that dEi n fij is smooth, since if we know 
(A.4) for smooth sets, then the general case follows by a rather standard approximation 
argument (which the weights do not greatly complicate). It is convenient to write 

E+ = Ei, E- = \ Ei, Ef = F-\Ef ) C E-^^i) =: a. 


Note that dEi is smooth in Ilj away from cAIj x {0}, and that dEi = E~^{dEi). Also, 
Hi has smooth boundary, so the relative isoperimetric inequality guarantees that there is 
some C = C{^li) < oo such that 


min?^"+™ 

± 


(Ef) < C 


n 


n+m—1 


{dE, n a) 


n-\-m 
n+m— 1 


(A.7) 


We will show below that, if 5 is small enough, then 

JE <2 in , (A.8) 

where JE denotes the Jacobian of F; and that, writing Eq := E\q^ 

^ < Jn-iFa in dE, (A.9) 

where Ju-iEq denotes the (n — l)-dimensional Jacobian, whose relevant properties are 
recalled below. We will first use these two facts to give the proof of (A.4), and then we 
give the details of (A.8) and (A.9). 

With this notation, we use the coarea formula to compute 


'E7 


p™ min f p^{x)dx^^^ c{m) min f ^{x))dx 

^ Je^ =*= Je^ 


(A.8) 


= c(m) min I JF{x,y)dxdy < C minT-T''^'^ [Ef), 


Similarly, using (more sophisticated) the coarea formula from Federer [18], section 3.2.22, 


, ~ ~ (A.9) 

rr+'^-^{dEinVLs) < c 


/ jFddrr-^ = c p^{x)dn 

oEi^ns JdEinn 


n—1 


Combining these two inequalities with (A.7), we deduce that (A.4) holds. 

Step 5. To complete the proof of the proposition, we verify (A.8) and (A.9). We fix 
{x,y) G and write M := DF{x,y). Then a computation shows that M = Mq + Mi 
where Mq, Mi are n x (n + m) matrices defined (in block form) by 


Mq = {Idn + (^ - l)V/9 (g) Vp, Vp (g) ^), 
P P 


Mi = (p-p)(v2p,0). 


This means that for V = (v,w} G M” x ME (understood as a column vector). 


MqV = V + 


(^-l)(u-Vp) + ^ 

P P 


Vp. 


(A.IO) 


One easily checks from the definitions that that |p — p| < J in so 

\MiV\<5\\V^p\\^\v\=C5\V\, 


(A.ll) 
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It also follows directly that 


kerM q = {(ctV p{x),w) : ap{x) + w ■ y = 0}, 
and can also check from (A. 10) that 

Mq is a linear isometry of (kerMo)"*" onto (A-12) 


Indeed, if we fix an orthonormal basis {fi,..., Vn} of M” such that vi = V/9(x), then the 
vectors dehned by 

Vi = —^—(p(x)Vp(x),y), Vi = {vi,0) for i = 2,... ,n, 

P[x,y) 


form an orthonormal basis for (kerMo)"*" that satisfies MoV) = Vi for all i. We now let 
14 + 1 , • • •, I4+m be an orthonormal basis for ker Mq, so that is an orthonormal 

basis for The Jacobian is 


1/2 


JF{x,y)=\ Y. \MVi,^...^MVi^ 


1/2 


(A.ll) 


^ I Mo 4*, A... A Mo 4*, 

l<2i 


+ 0{5) =l + 0{5). 


Thus (A.8) follows by a suitable choice of S. To verify (A.9), we compute Jn-iFg(x, y) 
for some aribitrary (x, y) £ dE. By dehnition, 


Jn—lFQ{x, y) 


E 


IMVi, A ... A MFi. 




1/2 


where ^ form an orthonormal basis for T(^^y^dE. which must contain 

kerMo(x,i/), since this is just the tangent space at (x,y) to E~^{F{x,y)) C dE. Thus 
dE has the form 

T(^x,y)dE = ker Mo © L, for some n — 1-dim. subspace L C kerM 5 '', 

And according to (A.12), Mq is a linear isometry of L onto its image. So we can choose 
a basis such that Mq!^ = 0 for i > n and {MqV}^Y form an orthonormal 

set. Then as above we find that Jn-iFQ{x,y) = 1 + 0{5), and we conclude by a suitable 
choice of <5. 

□ 
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